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Abstract 

The effect of spanwise-periodic mean-flow distortions (i.e. streamwise- vortex structures) 
on the evolution of small-amplitude, single-frequency instability waves in an otherwise 
two-dimensional shear flow is investigated. The streamwise-vortex structures are taken 
to be just weak enough so that the spatially growing instability waves behave (locally) 
like linear perturbations about a uni-directional transversely sheared mean flow. Nu- 
merical solutions are computed and discussed for both the mean flow and the instability 
waves. The influence of the streamwise-vortex wavelength on the properties of the most 
rapidly growing instability wave is also discussed. 

1. Introduction 

The effect streamwise-vortex structures have on the development of otherwise two- 
dimensional shear flows is of much interest for both theoretical and technological reasons. 
The theoretical interest stems from the role steady streamwise vortices play in amplify- 
ing three-dimensional disturbances in laminar shear flows. Transition to turbulence in rel- 
atively two-dimensional shear flows is virtually always preceded by the rapid growth of 
three-dimensional disturbances. The technological interest arises from the ability of steady 
streamwise-vortex structures to significantly enhance the mixing process. Mixing enhance- 
ment in shear flows has received renewed attention in recent years primarily due to the 
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current interest in supersonic combustion and noise reduction for the high-speed civil trans- 
port. The efficiency of scramjet propulsion systems depends on the time required to achieve 
a flammable (near-stoichiometric) mix of fuel and air within the combustor. Noise reduction 
in the mixer/ejector nozzle concept is accomplished by mixing entrained ambient air with 
high-speed exhaust plumes inside an acoustically treated ejector. 

Steady streamwise-vortex structures can arise in laminar-turbulent transition experi- 
ments due to wall roughness, free-stream disturbances and/or imperfections in the relatively 
two-dimensional excitation devices. In technological applications, the vortex structures may 
be produced by surface geometry or three-dimensional devices such as ramp injectors, tabs, 
or trip wires. In either case, the resulting three-dimensional mean flow will usually have a 
primary-flow direction - provided, of course, the cross flow associated with the streamwise- 
vortex structure is not too strong. Since the mean velocity component in the primary-flow 
direction is much larger than the velocity components in the cross-flow plane, the mean flow 
is governed by a set of equations that become of nearly boundary-layer type in the high- 
Reynolds-number limit (Davis k Rubin 1980). Small-amplitude disturbances propagating 
on such flows will behave (at least locally) like linear perturbations about a uni-directional 
transversely sheared mean flow. This paper is an attempt to provide a systematic first prin- 
ciples analysis of the linear stability characteristics of such disturbances using a combination 
of numerical techniques and high-Reynolds-number asymptotic methods. 

It is well known that steady streamwise-vortex structures can strongly influence the de- 
velopment of otherwise two-dimensional shear flows by amplifying three-dimensional distur- 
bances through a kind of parametric-resonance or secondary-instability mechanism (Nayfeh 
1981; Henningson 1987; Bennett k Hall 1988; Nayfeh k Al-Maaitah 1988; Hall k Seddougui 


2 



1989; and others). This phenomenon has most often been investigated in the context of the 
secondary instability of Gortler vortices on curved-wall boundary layers (e.g. Hall &; Horse- 
man 1991 and Li & Malik 1995). These investigations have been primarily concerned with 
determining the temporal amplification rates of the so-called fundamental secondary instabil- 
ities, i.e. the instabilities having a spanwise wavelength equal to that of the Gortler vortices. 
However, Li & Malik (1995) also investigated instabilities having spanwise wavelengths equal 
to twice that of the Gortler vortices and found that the temporal growth rates of these so- 
called subharmonic secondary instabilities can be comparable with the growth rates of the 
fundamental instabilities. 

Goldstein & Wundrow (1995, hereafter referred to as GW) analyzed the effect of intro- 
ducing a small-amplitude spanwise-periodic cross flow into a Blasius boundary layer. The 
imposed cross flow produces a distortion in the streamwise velocity component that initially 
grows with increasing downstream distance, reaches a maximum and then eventually decays 
through the action of viscosity. This decay occurs slowly and allows the distortion to desta- 
bilize the flow to three-dimensional inviscid disturbances over a relatively long streamwise 
distance. Goldstein & Wundrow concentrated on the limit where the spanwise wavelength 
of the mean-flow distortion is large compared to the boundary-layer thickness. This long- 
wavelength limit allows analytic solutions to be obtained for both the steady and unsteady 
components of the flow and therefore elucidates the relevant physical mechanisms. The au- 
thors also solved the spatial stability problem rather than the temporal one since it is the 
spatial theory that correctly describes convectively unstable flows. They found that the most 
amplified instability has a spanwise wavelength equal to twice that of the mean-flow distor- 
tion (i.e. it is a subharmonic secondary instability) and, furthermore, that the fundamental 
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secondary instability remains damped in the long-wavelength limit. 

The primary objective of the present paper is to extend the linear stability analysis 
of GW to the case where the spanwise wavelength of the mean-flow distortion is of the 
same order of magnitude as the shear-layer thickness. The linear stability analysis will also 
be generalized to spatially growing disturbances having an arbitrary, bounded behavior in 
the spanwise direction. This will enable an investigation of the effect of streamwise- vortex 
wavelength on the spatial growth rate and spanwise wavelength of the most rapidly growing 
secondary instability. The results of such an investigation should be of much technological 
interest since they could help guide the development of more efficient mixing-enhancement 
techniques. In order to carry out the investigation, numerical schemes must be developed for 
computing the evolution of the steady three-dimensional base flow as well as for determining 
the linear stability characteristics of a general time-periodic perturbation. These numerical 
schemes will also permit an evaluation of the follow-on nonlinear stage in the development 
of the most amplified instability wave analyzed in Wundrow & Goldstein (1994). 

The paper is organized as follows. In §2, governing equations for the steady base flow 
and the linear unsteady perturbation are derived in the high-Reynolds-number limit. The 
numerical schemes used to solve these equations are described in §§3 and 4. In §5, a compari- 
son of the results predicted by the long-wavelength solution developed in GW (and extended 
in appendix B) to those determined by a numerical solution to the aforementioned governing 
equations is presented. Numerical results for an order-one spanwise wavelength solution are 
also presented and discussed in §5. 
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2. Formulation 


Attention is focused on an incompressible shear flow formed at the interface between 
two parallel streams of differing velocity or, alternatively, between a single parallel stream 
and a flat plate. The Cartesian coordinate system ( x , y, z ) is attached to the interface with 
x in the direction of the external flow, y normal to the interface, and z in the spanwise 
direction. The origin of the coordinate system is located a distance L, downstream of the 
initial point of contact between the upper and lower streams or between the upper stream 
and the plate. All lengths are non-dimensionalized by 6, where 8, characterizes the local 
shear-layer thickness at x = 0. The * subscript is used to indicate dimensional quantities. 
The time t, velocity u = iu + jv + kw, and pressure p are non-dimensionalized by 8,/U., 
U, and pJJ}, respectively, where U, characterizes the velocity of the external flow and p, is 
the density. With this non-dimensionalization, the Navier-Stokes equations become 


it( + u- Vu + Vp = fl 1 V 2 u, 

(2.1) 

V* u = 0, 

(2.2) 

where V = id/dx + jd/dy + kd/dz is the gradient operator, 


R = hJJ.lv, > 1 

(2.3) 

is the local Reynolds number, v, is the kinematic viscosity and an independent variable used 
as a subscript denotes differentiation with respect to that variable. 

The solutions to (2.1) and (2.2) of interest here are expressed as the sum of a steady 
base flow and a time-dependent perturbation, 

U = U(x) + £u(x, t ), 

(2.4) 
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p = P(x) + ep{x,t) 


(2.5) 


where t < 1 characterizes the local amplitude of the perturbation at x = 0. The steady, 
spanwise-periodic, base flow {U, P} evolves over the long streamwise scale, 

x 2 = x/R, (2.6) 


and expands like 

U = iU 0 {x 2 , y,z) + R~ 1 V 0 (x 2 ,y,z) H , (2.7) 

P = P e {x 2 ) + • • ■ + R 2 Po(x 2 , y, z) + • • - , (2-8) 

for order one x 2 , y and 2 as R -> 00 , where V denotes the base-flow velocity in the transverse 
(or y-z) plane, P e denotes the pressure imposed on the shear layer by the external inviscid 
flow and Po is the first term in the pressure expansion to have a transverse variation. The 
unsteady flow {u ,p] also evolves on the x 2 scale, starting as a linear perturbation and 
becoming nonlinear sufficiently far downstream (for details see Wundrow & Goldstein 1994) . 
However, the present analysis only considers the local linear solution (i.e. the linear solution 
valid over streamwise distances of order <5„) so {u, p] are expressed as 


u = Re Jtt( 2 /, 2 )e 1 '*] 4- • • 

(2.9) 

P = Re [p(y, z)e‘*] + • ■ • , 

(2.10) 

for order one x, y, z and t, where 


X = R f a{x 2 )dx 2 - St, 

(2.11) 


a is the local streamwise wavenumber, S = 8 m F m /U* is the local Strouhal number (or non- 
dimensional angular frequency) and the dots indicate higher-order terms in e as well as 
R- 1 . 
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It was shown in GW that the leading-order base-flow solution is determined by the 


so-called boundary-region equations, 


U 0 U 0x2 + Vo-VrUo + K = V^C/o, (2.12) 

UoV 0x2 + Vo-V t V 0 + V t Po = V 2 t Vo, (2.13) 

U 0X2 + V T .V 0 = 0, (2.14) 


(Davis & Rubin 1980) where V r = jd/dy + kd/dz is the gradient operator in the transverse 
plane and a prime denotes differentiation with respect to the argument. The pressure fluc- 
tuation associated with the local instability wave was shown in GW to satisfy a generalized 
Rayleigh stability equation 


V T - 


' Vt P 

.(%-c ) 2 


° 2 P . =0 

(Uo - c) 2 U ’ 


(2.15) 


where c = S/a is the local phase speed. 

Equations (2.12)-(2.15) must be solved subject to 


{Uo, V 0 , p y } = 0 
{Uo, Vo, p} {C/e, jV 0 (+) , 0} 
for wall-bounded flows, or 


at y = 0 
as i/4 +oo 


(2.16) 


{U 0 , V 0 , p} {(C/ e 2 ± 1) 1/2 , jV 0 (±) , 0} as ±oo, (2.17) 


for unbounded flows, where U e {x 2 ) satisfies 

U e Ue' = -Pe', (2-18) 

and E 0 (+) and denote the scaled normal velocity components induced by the displace- 
ment thickness at the upper and lower edges of the shear layer, respectively. It will be 
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assumed that the spanwise varying part of the base flow remains confined to the shear layer. 
The external inviscid flow then remains two-dimensional to the required level of approxima- 
tion and the functions Vq ± \x 2 , y) are determined as part of the solution to (2.12)-(2.14). 

In addition to the normal-boundary conditions given by (2.16) or (2.17), (2.12)-(2.15) 
must be solved subject to boundary conditions in the spanwise direction. The steady base 
flow is required to be periodic in the spanwise direction with period 2n/{3. The spanwise- 
boundary condition for the unsteady perturbation is somewhat more general in that p need 
only remain bounded as z — > ±oo. 

3. Steady base flow 

The steady base flow is computed with a method similar to that used by Rubin, Khosla 
& Saari (1977) to analyze laminar flow in rectangular channels. The transversely varying 
component of the pressure is eliminated from the problem by applying the operator V T x 
to (2.13). This leads to 

Wo** + V 0 .V T U 0 + K = V* U 0 , (3-1) 

UqQq X2 + Vo-Vt-I^o — U 0x7 Qo — Vox 2 Uoz + W 0x2 u 0y — o, (3.2) 

00 = W 0y ~ Vo z , (3.3) 

u 0x7 + V T .V 0 = o, (3.4) 

where Qq is the leading-order term in the large-/? expansion of the streamwise component 
of vorticity. In view of (2.16), (2.17) and the requirement that the spanwise varying part of 
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the base flow remain confined to the shear layer, (3.1)-(3.4) must be solved subject to 

{U 0 , Vo, W 0 , A) - W 0y } = 0 at y = 0 
{Uo, Vo, Wo y , Qq} {U e , V 0 (+) , 0, 0} as y — > +co 
for wall-bounded flows, 

{Uo, Vo , VF 0y , R 0 ) -J- {(^e 2 ± 1) 1/2 , ^o (±) . 0> 0} as y -+ ±oo, (3.6) 

for unbounded flows, and, in either case, the solution must be periodic in 2 with period 

2n/p. 

In order to provide sufficient numerical resolution, the physical plane ( x 2 , y, z) is mapped 
into a computational plane (£, 77, Q using the coordinate transformations 

x 2 = Z-Z°, V = + sr = /3 _1 C (3-7) 

where £° = L*/R6 » is the scaled distance to the virtual origin of the shear flow, 0 = \Z£JU~ e 
accounts for viscous spreading, g is used to concentrate grid points near the dividing stream 
surface and 7/? allows for an arbitrary displacement of that surface. The specific forms of g 
and ip used in the computations are given in §5. For now, it suffices to note that y(0) = 0 and 
y(±oo) = ±00. Substituting (3.7) into (3.1)-(3.3) and introducing the rescaled dependent 
variables 

U 0 = U e U, V 0 = 0~ l V + y( ,Uo, W 0 = (3W, R 0 = R - y(U Qz , (3.8) 

yields, after some manipulation, 

DU = 2 '£UUz + - /r, (3.9) 

DR = iUR i - + ^Ur,W ( - 2 {K d g + k^UU^ + l -^U n c , (3.10) 

g y 
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(3.11) 


n = -w v - v c , 

9 

Wt + ±±£u + -,v v + 0 2 w< = o, 

^ 9 


(3.12) 


where 


(3.13) 


Ds ^(^- V ) + ^(|- H ')' 

» = 0 2 U' e is the pressure-gradient parameter, f3 = 00 is the spanwise-wavenumber parame- 
ter, and k$ = 9 3 Ug6" and = 6 3 U 3 i}>" account for curvature effects associated with the 
displacement thickness and dividing stream surface, respectively. 

Since (3.9)-(3.12) are parabolic in £, the base flow can be computed with a marching 
procedure starting at the streamwise position £ = £°. The solution is advanced in £ using a 
combination of interpolation and finite-difference formulas, 


f i+e = f(? + eA£, T 7 , C) = e/ i+1 + (1 - e)f + 0( A£ 2 ), (3.14) 

/* +£ = Me + ^ V , 0 = Ar'[(f + |)/ i+1 - 2 ef + (£- ^r 1 ] + 0( A£ 2 ), (3.15) 

which reproduce most of the familiar second-order accurate explicit and implicit schemes 
with appropriate choice of the parameter e. In particular, a semi-implicit Crank-Nicolson 
scheme is obtained by setting £ = 5 and a strongly implicit three-point-backward scheme 
results when e = 1. All the computations in the present investigation were done with e = 1. 

The coupled set of nonlinear equations that determine {U t+1 , V ,+1 , W ,+ 1 , 12 t+1 } are 
linearized using a Newton-Raphson procedure, 


Ui+e,n+l _ + ^ (3.16) 

f/ !+£ ’ n+1 _ yi+ £ ,n + + $ (3.17) 

9' 
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(3.18) 


w i+£ ' n+1 = w ,+€ ’ n + 8^- 

fli+e,n+l _ fli+s, n + ^ ( 3 . 19 ) 

where the superscript n indicates the iteration number and, following Rubin et al. (1977), 
corrections to the velocity components in the transverse plane are expressed in terms of a 
velocity potential 8$ and a stream function 8$. Introducing (3.14)-(3.19) into (3.9)-(3.12) 
and linearizing about {8U, 8^,8W, £12} = 0 gives 


[d* +£ - (^ £+ V +£ + 1 -1- Zfi i+S ^j U i+e 2 £ ,+£ t/' +£ ] 8U = FU i+£ , 

(3.20) 

^d* +£ 2s + l z i+£ u i+s ^ sn = f& +£ , 

(3.21) 

A£ (L i+£ <S!F + <S12) =Fr +£ , 

(3.22) 

A£L i+£ 54> = F # ,+£ , 

(3.23) 


(3.24) 


and the forcing functions FU, FQ , F$ and F# are given in appendix A. The components 
of the transverse velocity, V l+£ and VP +£ , in (3.20) and (3.21) have been lagged in the 
iteration (i.e. the 8<P and 8$ terms have been omitted) so that these equations decouple 
from (3.22) and (3.23). This allows an updated result for U l+C to be computed from (3.20) 
prior to solving (3.21)-(3.23). Consequently, the 8U terms in (3.21) and (3.23) have been 
omitted. The 8Q term in (3.22) has been retained because solving (3.21) and (3.22) as a 
coupled system greatly improves the convergence of the Newton-Raphson procedure for wall- 
bounded flows (Briley k McDonald 1984). In view of these considerations, the superscript n 
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has been dropped with the understanding that the latest iterates are used to evaluate terms 
at £ = £ ,+£ • 


(3.25a) 


Substituting (3.16)-(3.19) into (3.8) and the result together with (3.7) into the boundary 
conditions (3.5) and (3.6) shows that (3.20)-(3.23) must be solved subject to 

{U, v, W, g'Q -W v }= 0 at t) = 0 

{U, V , W v , Q) -4 {1, 0, 0} as 77-4+00 

and 

{6U, 69 v , g'89 c -6%, 69}= 0 at 77 = 0 
{6U, 8$r,(, 89, 8Q) -4 0 as 77 -4 +00 

for wall-bounded flows, 


(3.256) 


{U, V, W v , Q) -+ {(1 ± t/ e ' 2 ) 1/2 , F (± \ 0, 0} as 77-4 ±00, 


(3.26a) 


and 

{8U, 89 v ^, 69, 8fi} -4 0 as 77-4+00, (3.266) 

for unbounded flows, and, in either case, {U,V,W,Q} and {8U, 89, 89, 8Q} must be peri- 
odic in C with period 2 tt. The functions 77) appearing in these equations denote the 

rescaled normal-velocity components induced by the displacement thickness and are deter- 
mined as part of the solution to (3.20)-(3.23) (see appendix A). 

In view of the spanwise periodicity of the mean flow, the dependent variables in (3.20)- 

(3.23) are approximated by truncated Fourier series of the form 

+M 

/(fU c)«! £ (3-27) 

m— — M 

where 

fm = - r"/e- im< dC- (3.28) 

ft J — n 
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Derivatives in the normal direction are approximated with second-order-accurate central 
differences so, 

+ At? 2 ), (3.29) 

fvvm,j = = Ar? _2 (/^ j+1 - 2f' mJ + +0 {At) 2 ), (3.30) 

where Tfj = jArj. By introducing these approximations into (3.20)-(3.23), one obtains a 


system of algebraic equations that are expressed in matrix form as 

BDjSU + DU jSUj + ADjSU j+1 = FUf £ , (3.31) 

BDjSHj- i + DGjSiij + ADj6fij+i = F12’ +£ , (3.32) 

BLjS 1 + DLjS 9j + ALjS t j+1 -I- = FlP’ +£ , (3.33) 

BLjSZj - i + DLjStj + ALjS$j + i = F<P j +£ , (3.34) 


where <Jt/j = {6£7 m j}, FU j = {FU m j } , . . . are column vectors of length iV — 2AF -I- 1 and 
the non-zero elements of the order-N band matrices, BD V DU j , DQj, ADj, BLj, DL : 
and ALj, are given in appendix A. 

The LAPACK subroutine CGBSV is used to solve (3.31)-(3.34) subject to the boundary 
conditions obtained by introducing (3.27) and (3.29) into (3.256) or (3.266) (see appendix A). 
As mentioned above, (3.31) is solved first so that an updated result for can be computed 
from (3.16) and used in the remaining equations. Equation (3.34) is solved next. For wall- 
bounded flows, this provides the result for S $ needed to evaluate the 5 9 boundary conditions 
and so enables the solution of the coupled system (3.32) and (3.33). For unbounded flows, 
the 59 boundary conditions do not couple (3.32) to (3.33) so these equations can be solved 
seriatum. Updated results for V*J, and are then computed from (3-17)— (3-19) 
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and the solution procedure is repeated until the corrections {6U, <5<P, 8Q} become less 

than a preset tolerance. 

4. Unsteady perturbation 

The rescaled equation that governs the shape function p is obtained by substituting the 
coordinate transformation (3.7) together with the rescaled variables (3.8) into (2.15) to get 

(U -c) -If + P 2 p c< - a 2 p - UnPr, - 2p 2 U c p c = 0, (4.1) 

_9 Sri J 9 

where a = dot and c = cfU e are the rescaled streamwise wavenumber and phase speed, 
respectively. In view of (2.16) and (2.17), (4.1) must be solved subject to 

p v = 0 at 77 = 0, p 0 as 77-* +00, (4.2) 

for wall-bounded flows, or 

p -¥ 0 as 77 — ► ±00, (4.3) 

for unbounded flows. 

Since interest here is in disturbances that remain bounded as z ±cc and since the 
mean flow is periodic in the spanwise direction, Floquet theory can be used to express the 
relevant solutions to (4.1) as a summation of solutions of the form 

P = p(v, Cl7)e l5< , (4.4) 

where p is periodic in £ with period 2 7r and 7 is a real characteristic exponent. Substituting 
(4.4) into (4.1), approximating U and p by truncated Fourier series of the form (3.27) and 
equating the coefficients of like powers of e l( » yields 

min (M,m+M) . T 1 / 1 \ 

Y. ] (Um-n ~ 2 c5 m , n ) — {'"/PnrjJ ~ ^ nPn 

n=max(— M,m— M ) ^ 
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for — M ^ m ^ +M, where 


- -^U m -nr,Pnr, + 2/? 2 (t7 + 7 ) (m - n)[/ m _ n p n j = 0 , (4.5) 


\ m = [a 2 + /? 2 (m + 7) 2 ] ^ , ReA m > 0, (4.6) 

and <5 m>n denotes the Kronecker delta. It follows from (3.25a), (3.26a), (4.2) and (4.3) that 
(4.5) must be solved subject to 

p m „ = 0 at 77 = 0 , p mv + \ m g'pm^Q as 77-M-oo, (4.7) 


for wall-bounded flows, or 

Pm v ± Xmg'pm -> 0 as 77 -> ±00, (4.8) 

for unbounded flows, where these expressions are obtained by using the exponential behavior 
of p m in the free stream. 

Introducing the central-difference approximations (3.29) and (3.30) into (4.5) leads to 
a system of algebraic equations that can be written in matrix form as 

bjpj-i + djpj + ajp j+1 = 0, (4.9) 


where pj = {p mj } is a column vector of length N = 2 M + 1 and bj, dj and aj are band 
matrices of order N. Away from the boundaries, the non-zero elements of bj, dj and a ; are 
given by 


1 1 


- ~£~2 ^ + 2 gi. j gi 2 (Um-nj ^m,n) + ^ ^, 2 U m -n n ■ , ( 4 ’ 10 ) 


2 1 


d m , n J = - ( + x nj ~ 2c6 m , n ) + 2f(n + j)(m-n)U m -n,j, (4.11) 
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2 1 


(4.12) 


n,n,j — ^^2 gi 2 (Um-n,j 2c$m,n) bm,n,ji 

where -M ^ m ^ +M and ma x(-M, m - Af) ^ n ^ min(M,m + M). Expressions for 
b m ,n,j, ^m,n,j and a m>nj at the boundaries are obtained from (4.7) or (4.8) by using (3.29) 
to approximate the r) derivatives. 

For a given base-flow profile U, Strouhal number S and characteristic exponent 7 , the 
rescaled wavenumber a is determined by solving the eigenvalue problem (4.9). This is done 
with an iterative method similar to that used by Hughes (1972) to solve the Orr-Sommerfeld 
equation. The method is most easily applied when (4.9) is rewritten as 


A x — 0, 


(4.13) 


where = P m ,j an< ^ the non-zero elements of A are given by 

A m +Nj,n+Nj—N = bm,n,ji A m +N j,n+N j n+Nj+N ( 4 - 14 ) 

A normalization condition, Xk = 1 at k = k+, is imposed on x and the A:*-th equation in 
(4.13) is relaxed by allowing 

5>*..*** = d^ 0. (4.15) 

k 

Equation (4.13) is then recast as an inhomogeneous problem by moving the terms involving 
a:** to the right-hand side and replacing the fc*-th equation with the identity = 1. This 
yields 

Bx = /, (4.16) 

where 

Bk, i = + (1 - <***,/)(! - fk = &k, k* - (1 - <5fe,fc*)Afc,^. (4.17) 
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Given an estimate for a, x is determined by solving the sparse linear system (4.16) with the 
LAPACK subroutine CGBSV. The eigenvalue d is then found by locating the zeros of A 
using Newton’s method. 

5. Results and discussion 

The generalized Rayleigh stability problem defined by (4.1)-(4.4) or equivalently (4.1)- 
(4.3) and the spanwise-boundary condition, 

p(v > C + 2 7r 1 7) = e l2K7 p(77, C I t)i (5- 1 ) 

has a number symmetries with respect to the characteristic exponent 7 that are used to sim- 
plify the numerical computations. Let {p(rj, C 1 7), 0(7)} denote an eigenfunction/eigenvalue 
pair determined for a particular value of 7, it then follows from (5.1) that there exists an 
infinite sequence of eigenpairs {p(rj, Cl 7 + k ) , a(j + k ) } such that 

p{v, C I 7 + ^) = PiV’ C I 7) and a(7 + fc) = 0(7), (5.2) 

where k = 0, ±1, ±2, If the base-flow profile U is an even function of £ as it will now be 

assumed, there also exists an eigenpair {^(77, (| - 7),a(-7)} suc h that 

p{r),-0 -7) = ?(»7»Cl7) and a(-7) = a(7). (5.3) 

Combining this result with (5.2) leads to 

p{T}, — C 1 1 — 7)=P0?,Cl7) and a(l - 7) = 0(7). (5.4) 

Therefore, only solutions for 7 in the range 0 to \ need be computed since solutions outside 
this range can be easily obtained from (5.2)-(5.4). 
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The C symmetries of the solutions to the generalized Rayleigh stability problem are also 


of interest. Let the eigenfunction p(rj, £ 1 7) be expressed as 


Piv, Cl 7) =Po(*7,Cl7) + Pe(*?,Cl7), ( 5 - 5 ) 

where p 0 and p e are odd and even functions of (, respectively. Then, since U is an even func- 
tion of C) Po and pe satisfy (4.1)-(4.3) independently but are coupled through the spanwise- 
boundary condition (5.1). This coupling becomes clear when (5.5) is used to rewrite (5.1) 
as 

Po{ 7?,C + 2tt |7) = cos(2Tr7)p 0 (7 ? ,C|7) + isin(27t7)p e (77,Cl7), (5-6a) 


and 


Pe{vX + 2 tt I 7) = i sin (2 717 )^( 77 , C I 7) + cos(2 ^^)p e (j}^ C I 7)* (5.66) 


It follows that the coupling terms in (5.6) drop out whenever 7 = \k. p 0 and p e are then 
independent solutions to the generalized Rayleigh stability problem which (together with 
their corresponding eigenvalues) are termed odd- and even-mode solutions. 

The numerical computations discussed below are for wall-bounded flows so the func- 
tion ^(f) introduced in (3.7) is identically zero. The stretching function g(rj) used in the 
computations is given as 


/ArcsinhAr? . _ 

g = ■ + 1 B 

\ sinh AAtj 


tanh (c 


A 77 sinh Arj\ 
sinh A At] J ’ 


(5.7) 


where A, B , and C are real constants. Typical values for these constants are A ^ 5.5 x 
10“ 6 A?7^ 1 , B = -0.4 and C — 10 when the computation is done off the real axis and 
A ^ 5.5 X 10" 6 At/“ 1 , 5 = 0 and C = 00 when it is done on the real axis. 

Figure l is a comparison of the scaled growth rate —6a{/a A predicated by the long- 
wavelength solution developed in GW (and extended in appendix B) to that determined 
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from a corresponding numerical solution to the generalized Rayleigh stability problem for 
the analytic base-flow profile 

U = U B (g) - a 3 ge- g3 ^ 3 cosC, (5.8) 

where 0 < a < 1 is a scale factor and U B is the Blasius profile determined by 

U B - V B U' B = 0 and \U b + V' b = 0, (5.9) 

with U B ( 0) = Vs(0) = 0 and U B {oo) = 1. The computations were performed with the grid 
parameters At? = 0.001, J = 1500 and N = 27. The curves are plotted for both odd- and 
even-mode solutions at y = i. At the largest value of cr, the exact and asymptotic results 
appear quite different. The exact solutions show a relatively fast growing odd mode and 
a weakly unstable even mode whereas the asymptotic results predict a completely damped 
odd mode and a rapidly growing even mode. However, this disparity diminishes and the 
exact solutions approach their respective asymptotes as a — > 0. 

Figure 2 is a plot of the scaled growth rate —0ai/a 4 as a function of the characteristic 
exponent 7 for the analytic base-flow profile (5.8). The curves are plotted for a scaled 
Strouhal number 9S/cr 2 U e = 0.06 which is near the peak in the asymptotic growth rate of 
the even mode at 7 = i. Also shown are the growth rates predicted by the long- wavelength 
solutions (B46) and (B47) which have been evaluated with a — 0.1. The results plotted 
in parts a and b of the figure correspond to solutions that become, respectively, odd and 
even functions of C when 7 = As noted above, the odd and even mode designations 
apply only when 7 = \k. Indeed, these designations change as 7 is decreased from \ to 0, 
i.e. the solution that corresponds to an odd mode at 7 = | becomes an even function of £ 
when 7 = 0 and vice versa for the solution that corresponds to an even mode at 7 = 


20 





Figure 2: Scaled growth rate vs. characteristic exponent for U given by (5.8), 6/3/a = 
6S/a 2 U e = 0.06 and various values of a. Solid lines, exact numerical solution; dashed lines, 
long-wavelength asymptotic solution, (a) solution that is odd mode at 7 = (5) solution 

that is even mode at 7 = |. 
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As in figure 1, the exact and asymptotic solutions in figure 2 appear quite different at the 
largest value of a but come into agreement as a decreases. It is interesting to note that the 
solution for a = 0.2 shown in figure 2 a is damped at 7 = 0 and | but is unstable when 
7 is near 0.35. This points out the need for investigating more than just the fundamental 
and subharmonic solutions to the generalized Rayleigh stability problem. Figure 2 also 
shows that the effect of the streamwise-vortex structure on the instability waves becomes 
concentrated near 7 = | as a becomes small which is consistent with the conclusions in GW. 
This means that, once a becomes sufficiently small, the base flow is inviscidly unstable only 
to subharmonic disturbances. 

Figures 1 and 2 are plotted for a scaled span wise wavenumber 60 /cr = i. In order to 
demonstrate the effect of this parameter on the linear-stability results, the neutral (or zero- 
growth) values of the scaled Strouhal number 6S/a 2 U e predicted by the long-wavelength 
solution are plotted against 60/a in figure 3 for both odd- and even-mode solutions at 
7 = i. The figure shows that there are multiple ranges of 60 /a over which the profile (5.8) 
is inviscidly unstable. However, each range is finite so above a certain value of 60 /a (5.8) 
becomes inviscidly stable even though it remains inflexional over a portion of the spanwise 
domain. This is consistent with the experimental observations of Hamilton & Abernathy 
(1994) who investigated the conditions under which transition to turbulence can be caused 
by introducing streamwise vortices into an otherwise two-dimensional, non-transitional shear 
flow. Figure 3 suggests that growing solutions exist when 60/a — 0. This result must be 
disregarded because the analysis of GW shows that the long-wavelength asymptotic solution 
breaks down as 0 -4 0 with a > 0 held constant. In order to show that the above conclusions 
apply to other base-flow profiles, a plot of the neutral Strouhal number versus the spanwise 
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8S/a 2 U, 


9S/a 2 U. 




Figure 3: Scaled neutral Strouhal number vs. scaled spanwise wavenumber for the long- 
wavelength asymptotic solution with U given by (5.8) and 7 = 5. (a) odd mode; (6) even 


unstable 


unstable 


Figure 4: Scaled neutral Strouhal number vs. scaled spanwise wavenumber in the notation of 
GW for the long-wavelength asymptotic solution with 7 = \ and U given by (3.32), (3.43) 
and (3.51) of GW where a = 0.05, Bq= — 16 7r and A = 0. (a) odd mode; ( b ) even mode. 





wavenumber is shown in figure 4 for the base-flow profile derived in GW. 

The time needed to compute curves corresponding to figure 3 for the exact solutions 
to the generalized Rayleigh equation would be prohibitive. However, the results in figures 1 
and 2 suggest that, at small values of a, qualitatively similar curves would be obtained. 
It is unlikely that this similarity would persist as a increases but it is interesting to note 
that, in their numerical investigation of the secondary instability of Gortler vortices, Li & 
Malik (1995) found that the relative growth rates of the odd and even modes are strongly 
influenced by the Gortler- vortex wavelength (i.e. the spanwise wavelength of the base flow). 

The influence of non-parallel-flow effects is demonstrated by computing the local sta- 
bility characteristics of a base-flow profile U that satisfies the boundary-region equations 
(3.9)-(3.12). The computation is done for a wall-bounded flow with U e — 1 and 

U = U B , V = V B + 0 2 U B H o W = -U B Hg , n = -(U B H g ) g -p 2 U B H co (5.10) 

at f = 1, where U B and V B satisfy (5.9), P = y/£, and the function H(rj,Q is given as 

H = A^f/gt/gSinC, (5.11) 

with X B = U' B {0). The form of the ‘initial’ condition (5.10) has been chosen so that U, V , 
and W satisfy the equation obtained by eliminating the U £ terms between (3.9) and (3.12). 
Equations (3.9)— (3.12) are solved subject to (3.25) and (5.10) using the method described in 
§3 with = 0.001, At] = 0.005, J = 900 and N = 17. 

Figure 5 shows lines of constant streamwise velocity U in the plane at various 
streamwise stations. The contours at £ = 1 are not shown but, in view of (5.10), they 
would be straight lines parallel to the £ axis. The contours in parts a-e of figure 5 reveal 
a rapid departure from the initial two-dimensional state as the distortion in U produced 
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contours in the g-Q plane. Thick lines, U contours; thin lines. 
6) ? = 1.15; (c) £ = 1.2; (d) £ = 1.3; (e) £ = 1.5; (/) £ = 1.8 






Figure 5: Streamwise velocity contours in the <?-£ plane. Thick lines, U contours; thin lines, 
inflexion points, (a) £ = 1.1; ( b ) £ — 1.15; (c) £ = 1.2; (d) £ = 1.3; (e) £ = 1.5; (/) £ = 1.8. 
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by the imposed cross flow grows with increasing £. This type of disturbance growth is 
similar to the algebraic growth proposed by Ellingsen & Palm (1975), Hultgren & Gustavsson 
(1981) and Landahl (1990) as an alternative or ‘bypass’ transition mechanism. However, 
only the base-flow distortion undergoes this type of growth in the present study and the 
follow-on secondary instability (discussed below) exhibits the more conventional exponential 
growth. The wavy appearance of the contours in figure 5 is similar that observed in the early 
stages of Gortler- vortex development (Swearingen & Blackwelder 1987). The ‘mushroom- 
like’ structures characteristic of the later stages of that development do not appear in figure 5 
because the vortex structure shown there is not sustained by wall-curvature effects. The 
growth of the U distortion in the present study is eventually reversed by viscous-diffusion 
effects and the base flow tends towards its two-dimensional undisturbed state as £ becomes 
large (cf. figure 5/). 

Figure 5 also shows the location of the inflection points in the U-g plane, i.e. the points 
where U gg = 0. Wundrow & Goldstein (1994) showed that the connection between the in- 
flexions in U and the stability of the base flow is much more complicated for a transversely 
sheared base flow of the kind being considered here than for a plane flow for which Rayleigh’s 
inflexion-point theorem applies. Nevertheless the inflexions in U do indicate a redistribu- 
tion in momentum produced by the streamwise-vortex structure that does play a role in 
determining the local stability of the flow. 

Figure 6 is a plot of the local instability-wave growth rate — 0a, against the local Strouhal 
number 6S/U e for the base-flow profile U shown in figure 5. Growth rate curves are shown 
for odd-fundamental (7 = 0, odd-mode) and even-subharmonic (7 = |, even-mode) solu- 
tions. Curves for even-fundamental and odd-subharmonic solutions are not shown since only 
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damped forms of these solutions were found. It should be noted that the shooting method 
used here does not determine all the eigenvalues of the generalized Rayleigh stability problem 
and a global method must be applied to show that there are no unstable even-fundamental 
or odd-subharmonic solutions. In figure 6a, only the even-subharmonic solution is shown be- 
cause the odd-fundamental solution is damped at this streamwise station . This suggests that 
at short distances (on the z 2 scale) downstream of the point where the cross flow is introduced 
the linear instability of the flow is dominated by a subharmonic rather than a fundamental 
solution. This is also the case in the long- wavelength analysis considered in GW. Further 
downstream, as the U distortion becomes larger, the odd-fundamental solution becomes 
unstable and soon exhibits a peak growth rate larger than that of the even-subharmonic 
solution. The growth rates of both solutions continue to increase until £ « 1.3 after which 
point they decrease with increasing £. It is interesting to note that the streamwise position 
of peak growth precedes the position were the U distortion is a maximum. This behavior is 
also evident in the long- wavelength asymptotic solution developed by GW. 

Figures 5 and 6 show that the imposed cross flow (5.10) produces a streamwise region in 
which the boundary layer is inviscidly unstable. In order to show that this region of instability 
can significantly affect the boundary-layer development, an estimate of its streamwise extent 
is made using the flow parameters reported in Swearingen & Blackwelder (1987). In their 
experimental investigation of the evolution of Gortler vortices, Swearingen & Blackwelder 
introduced a small-amplitude acoustic disturbance 60 cm downstream of the leading edge 
of a concave plate and observed vortex breakdown leading to laminar-turbulent transition 
at downstream distances between 100 and 110 cm. Taking Z, = 60 cm and using £° = 
L./R5, = 1 together with (2.6) and (3.7), the dimensional distance from the leading edge 
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in the present investigation is 


x* + Z* = RS* £ = 60 £ cm. 

Therefore, since £ ranges between LI and 1.8 in figures 5 and 6, the distance over which the 
imposed cross flow (5.10) destabilizes the boundary layer is the same as the distance over 
which the secondary instability of Gortler vortices leads to laminar-turbulent transition in 
the Swearingen h Blackwelder (1987) experiment. 

The ability of a given streamwise- vortex structure to significantly affect the development 
of an otherwise two-dimensional shear flow will, of course, depend on the vortex strength. 
The strength of the streamwise-vortex structure can be characterized by its circulation, 

1 f n /& r°° 1 /*°° r 1 * 

rs rJ 0 L ( 5 - 12 ) 

where use has been made of (3.3), (3.5), (3.7) and (3.8). The non-dimensional circulation 
(or vortex Reynolds number) r of the imposed cross flow (5.10) is approximately 9 - well 
within the range of vortex strength considered in the experimental investigation of Hamilton 
fc Abernathy (1994). Vortex structures of lesser strength were also found to destabilize the 
boundary layer but the streamwise distance over which this occured tended to be smaller. 
Concentrating the imposed cross flow near the wall was found to reduce the streamwise 
distance needed to achieve an inflexional U profile but again the resulting vortex structures 
tended to damp out sooner. Changing the spanwise wavelength of the vortex structure was 
found to alter the linear stability characteristics however further investigation is need to 
verify if behavior such as that shown in figures 3 and 4 occurs in the order-one wavelength 
case. 
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The results presented here show that the introduction of even a very weak streamwise- 
vortex structure into an otherwise two-dimensional boundary layer can lead to the ampli- 
fication of certain three-dimensional disturbances through a kind of parametric-resonance 
or secondary-instability mechanism. The amplification rates of these secondary instabilities 
are strongly dependent on their spanwise wavelength as well as the strength and spanwise 
wavelength of the streamwise-vortex structure. Nevertheless, at sufficiently high Reynolds 
numbers, the secondary instabilities can grow at a rate much faster than that of the Tollmien- 
Schlichting waves which would otherwise dominate the initial stages of laminar-turbulent 
transition. Viscous effects eventually damp out the base-flow distortion and consequently 
any instability wave amplified by it. However, the damping process occurs over the relatively 
long X 2 scale so an initially linear instability wave can easily become nonlinear before this 
process is complete. 

Nonlinear effects will become important first within a thin critical layer located at the 
transverse position where the phase speed c of the instability wave equals the base-flow 
velocity U (once the instability-wave amplitude and growth rate become sufficiently large 
and small, respectively). The unsteady flow outside the critical layer remains essentially 
linear but the instability-wave amplitude is then completely determined by the nonlinear 
motion inside the critical layer. An analysis of this stage of evolution is given in GW for 
the long wavelength limit and in Wundrow & Goldstein (1994) for the order-one wavelength 
problem considered here. It is worth noting (as pointed out by Hultgren 1992) that the near- 
neutral approximation used in the order-one-wavelength critical-layer analysis will remain 
accurate even at frequencies close to the peak growth because the slope of the growth rate 
curve between that point and the upper branch is nearly constant (cf. figure 6). 
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Appendix A. Coefficients and boundary conditions for (3.31)— (3.34) 


In this appendix, expressions are given for the forcing terms, matrix elements and bound- 
ary conditions associated with the algebraic equations (3.31)-(3.34). The forcing terms are 
defined as 


FQ s -A£ 


FU = -A£ (DU - 2 ZUUs - i±^t / 2 + ^ , 

(Al) 

- tUfy + ZUtVf. - ^U n Wt + 2{Kgg + K^)UU C - —y-U V { 
9 9 

, (A 2) 

III 

1 

t> 

1 

| k- 

£ 

+ 

(A3) 

F$ = -AC (V* + l ~~^U + ±Vr, + , 

(A 4) 


where it should be noted that these terms vanish when {U, V, W , 12} satisfy (3.9)-(3.12). The 
non-zero matrix elements are given by 


BL m<n j — BD my n,j 




4 At? 5 ' 


r V m - n , 3 -l - 6m , 2 g ') g , 2 - 


DL n = DD m , nJ + f ?T +e i mWj£ nii = -«5 m ,„A^ (-1^ + ^rn^j , 


DUm,n,j — DD mn J 


Mi 1 yi+l 
4 A 7? p' 

& + i fl+ . + 1 + r - A J „ 


AL m,n,j - AD m ,n,j + AAr)n f V m-n,j + 1 “ Ar)2 ( 1 2 <?' j p'. 2 ’ 


A£ (, At?<?'^ 

At ? 2 


2 s 


Or! 1 

no — nn . _ Z f l+e f/ 1+£ . 


(A 5) 

(A 6 ) 

(A 7) 

(A 8 ) 
(A 9) 


where -M ^ m ^ +M, max(-M, m - M) ^ n < min(M, m-f M), S m , n denotes the 
Kronecker delta and dependent variables with superscript i + s are evaluated using (3.14) 
and (3.15). 
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Setting Atj = r} ma , x /J and introducing (3.27) and (3.29) into (3.25) and (3.26) shows 
that, in algebraic form, the normal-boundary conditions read 

Um , 0 = V"m,0 = Mm,0 = 0) Wjn,-1 = W m ,l ~ 

U m ,J = 2(5 m> 0, TTlV mi J = 0, W m ,J+i = ^m,J = 0) (A 10 fl ) 

and 

^ # m ,o = o, 5<p m - 1 = 6# m> i, = $*m,i - 2imA^# m ,o, 

5Um,j = 6<p m ,j = 612 mi j = 0, m5$ mi j+i - m5$ m ,J-u (A 106) 

for wall-bounded flows, and 

U m ,±j = 26 m ,o(l±t/ e -2 ) 1/2 , rnF m ,±j = 0, W m ,±j±i = W'm.iJTi , 12 m ,±j = 0, (A 11a) 

and 

&Um,±J = & $m,±J — $ ^m,±J = 0, ^m,± J^l i ( A 1 1 6) 

for unbounded flows. The normalization condition 

<5 #o,o = 0, (A 12) 

and, in the case of unbounded flows, the kinematic condition 

S$o,i = 6*o,-u (A 13) 

must also be imposed on the solution in order to uniquely determine <5# 0 j. These conditions 
allow the m = 0 component of (3.34) to be solved as an initial-value problem. With <5# 0 ,j 
so determined, the normal-velocity components induced by the displacement thickness 
can be found from (3.17). 
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By substituting (3.27), (3.29) and (3.30) into (3.18) and the ^derivative of (3.17) and 
combining the results with (A 10) and (All), the following supplemental expressions are 
obtained, 

V m ,.i = V mAi (A 14) 


6V m< J + 1 = 8V m> j - 1 + 1 - 6$m,j) - 1 - 6#m,J+ 1 ) 


-i roj5 2<+e M*w.J- 1 , 

(A 15) 

= ITTlS *P 7 n ^j -p 

(A 16) 

for wall-bounded flows, and 


bVm,±J±\ = 1 i ~ T /± (<^m,±JqFl 

- <5^m,±J±l) 


(A 17) 

= i g /± jh±j6 ^m,± ? 

(A 18) 

for unbounded flows, where 


8V=^-S^ n + 0 2i+e S^ o 

5 9 

4 4 g r ±J 

,±J = A V g- ±] ' md h±J -2g' ±J *A wlj 


These expressions are used when computing updated results for V^J and at the tj 


boundaries. 


Appendix B. Long- wavelength solution to (4.1) and (4.2) 

In this appendix, an asymptotic solution to the generalized Rayleigh stability problem 
(4.1) and (4.2) that applies in the limit as /3 ->■ 0 is constructed. The analysis closely parallels 
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that given in GW but is carried out here for shape functions p having the general Floquet 
form (4.4). The relevant wavenumber and phase speed scalings for the long-wavelength limit 
are 

ft = a/3, a = era = aa r + i<7 4 &;, c = ac = ac r + icr 4 c,, (B 1) 

where 0 < a <C 1 is a scale factor, fi is an order-one real constant and the real quantities ci r , 
<5j, c r and c,- have expansions of the form 

a r = a r0 + aa rl + •••, (B2) 


as a — )■ 0. 

As in GW, the present analysis is restricted to wall-bounded flows that behave like 
small perturbations about the Blasius solution. In the streamwise region of interest, the 
streamwise velocity of the base flow is then expressed as 

U = U b + o 4 U d , (B 3) 


where U B corresponds to the Blasius solution determined by (5.9) and U D is a spanwise- 
periodic distortion. Introducing (B 1) into (4.1) and (4.2) leads to 


Pg 

_(U - <7C) 2 . 


+ a 2 P 


PC 


[(U-crcyit 


2 -2 

- a a 


( U - ac ) 2 


= 0 , 


(B4) 


and 


Pg — 0 at <7 = 0, p — > 0 as g — ^ -|-oo, (B 5) 

where, for convenience, g is treated as an independent variable. 

Since much of the analysis in GW carries over directly, only the final results are presented 
here. In the inviscid wall layer described by g = a~ l g = 0(1), the base flow is given as 

U = a\ B g + a 4 (if — ±\ 2 B g 4 ) + ■ ■ ■ , (B 6) 
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where A B = £7^(0) and U{g, C; o) = U D {ag, C; = 0(1). The shape function p is given by 


p = a A + 0(a 5 ) 


(B 7 a) 


and 


El = E1 V 

( U-CFCr ) 2 A| 


C7 6 

Iddg - 9c ) J + A| 


•f + /ity* -</>)- 


ic.ff 2 


< 7?(<7 “ 9c) 2 


+ (B 76) 


where .4 is an arbitrary function of Q that has an expansion in a containing terms up to but 
not including order cr 4 , V = 0 (A^d/dC + A^) - a 2 A, the normal position of the critical 
level, g c = cr£r c , is determined by the condition 


U{g c X;<?) = °Cr, 


(B8) 


the ± superscript indicates different values for g ^ g c , a subscript c denotes evaluation at 
g = g c , $ is a real function of g and £ defined as 


~to Hs-s ,) 3 


■f denotes the Cauchy principal value, and 


ft — Uggc 4 ^b9c 


(B 9) 


(BIO) 


is the scaled normal derivative of the base-flow vorticity at the critical level. The constant 
— <£“ is set equal to in which corresponds to a logarithmic phase shift of n across the 
linear critical layer. 

In order to obtain solutions for p of the form (4.4) that are valid for all 7 , it is necessary 
to retain some intermediate terms (which were not needed in the GW analysis) in the main- 
boundary-layer and outer-layer solutions. In the main boundary layer described by g = 0(1), 
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the solution for p is given, to the required order of accuracy, by 

p = crA + a 2 j\u B - ac r ) 2 [(1 - a 2 JC)VB - a (I - a 2 I< C)CA] dg + • • • , (B 11) 

where 5 is a function of ( that has an expansion in a containing terms up to but not including 
order a 4 , C = p 2 d 2 /d( 2 - a 2 and 7, J and I< are the real functions of g determined by 


(U B - ac r ) 2 I' = 1, and [{U B - ac r ) 2 { J, K} 1 ]' = ( U B - ac r ) {1, /}, 


(B 12) 


with {/, J, J', K, K'} = 0 at g = 0. Matching (B 11) with the wall-layer solution (B 7) shows 
that 

B = X- 1 c~ 1 +a 3 u, (B 13) 

where u> is a complex function of Q defined as 
u> = X B 3 lim \$ + \X 2 B {^f + \g c g + g 2 c \nag) + inp-^c 2 ^ + - ic,-), (B 14) 

g-*oo L 

and use has been made of the small-cr expansion 


5 c = A- 1 c r -a 3 (A- 1 f/ c -^A- 3 ^) + ---, 

which is easily obtained from (B 6) and (B 8). 

In the outer layer described by g = ag = 0(1), the solution for p expands like 

p = ap + 0(cr 5 ), 

where the function p of g and ( is determined by the Helmholtz equation, 


Pgg + ~ Oi 2 p = 0, 


(B 15) 


(B 16) 


(B 17) 


which must be solved subject to 


■ = (1 + av 2 C — o^v^C 2 ) A, pg = (z^i C — g 2 vzC 2 )A-\- a 3 Vcji at g — 0 } (B 18a) 
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in order for (B 16) to match with the main-boundary-layer solution (B 11) and 

p — > 0 as g — ► 4-oo (B 186) 

in order to satisfy the free-stream boundary condition. The v n appearing in (B 18a) are real 
constants defined as 

XgC r I ^ 2 i ^ 3 i ^ 4 } — { 1 gXqC-tI oo? J ^ (7\sC r K. 00 i J qg j -^'oo}'? (B 19) 

where 

/oo = lim (J - O), 43(1- o-c r )“ 2 , (B 20) 

5 — >°c 

Joo = lim (J - \g 2 - = lim {J' - g) , (B21) 

g ~+ 00 z 5—^00 

Koo = ^( K - k 1 '^ 3 - 2 J oo9 2 - O). - KooS J - 7 °°2) , (B 22) 

and 

roo 

Loo = / [(U B - ac r ) 2 4 J - \g 2 - .O - *Joo]d<?. (B 23) 

Jo 

The relevant solutions to (B 17) and (B 18) have the form 

+ oo 

{p,A}=| £ {Pme- Xm ' 9 ,A m }e'l m+ ™, (B24) 

m=— oo 

where 7 is the real characteristic exponent. It follows from (B 14), (B 17), (B 18) and (B24) 
that 

A m = [d 2 + /3 2 (m + 7) 2 ] 1/2 , ReA m > 0, (B25) 

Prn = (1 - av 2 \ 2 m - a 3 v 4 \n)A m , (B26) 

and 

6m -Tm — 1 4" d m Am 4“ + l (B 27) 
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where 


bm = <j 3 | u; l(^m -2 ““ ^m-1 “ ^m)i (B 28) 

<^m = A m — (l/i + |a 3 CJo) — az/ 2 ^m — ” a (B 29) 

a m = + 4A^ +1 — A^ +2 )> (® 

Wm = - [ +n ue- im <d<;, (B 31) 

7T J— tc 


and, in order to fix ideas, it has been assumed that U D is proportional to cos Equations 
(B25)-(B27) together with the definition of the scaled Strouhal number, 

S = ac=~4 7r S, ImS = 0, (B 32) 

O’ f/g 

form the eigenvalue problem that determines a as a function of 5, /? and 7 . 

The small-cr behavior of a is determined by using the method of cyclic reduction to 
express (B27) as 

— b m —ib m d m +i A m _2 + {dm-ldm.dm+1 ~ a m—lbmdm+l ~ 

A m +2 = 0; (B 33) 

and then imposing the normalization condition Ao = 1 so that attention is restricted to 
instability waves having the bulk of their energy in the m = 0 mode. Equations (B28)- 
(B30) together with (B33) show that A± 2 = 0(a 3 ) and therefore that 

cLido^i — ct—ibodi — d-iaobi = 0(<r 9 ), (B 34) 

which, to the required level of approximation, is the characteristic equation for (B27). 

The solutions for a that are of most interest here are those corresponding to instability 
waves having spanwise wavelengths nearly twice that of the base flow since it is in this 
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parameter range that the waves exhibit their most rapid growth. Substituting 


■y — 2 "b (B 35) 

into (B25) and the result into (B28)-(B30) yields 

b m = [k 2 m ~ ?{m+ |)]+0(a 6 ), (B 36) 

dm = k m - (i/1 + ±cr 3 u>o)kl - av 2 k 3 m - <J 2 v z k^ - c 3 v 4 k 5 m 

- a 3 (2m+l)(^ - | k-'tfh + Oicr 4 ), (B37) 

Gm = -a 3 ^ [P m +0 2 (m+ i)] + 0(<j 6 ), (B 38) 

where 

k m =[a 2 + 0 2 (m + i) 2 ] 1/2 , Rek m >0. (B39) 


Then, since k - 1 = £ 0 and the m = 0 component of (B27) implies do = 0(cr 3 ), substitution 
of (B36)-(B38) into (B 34) leads to 

[jfco - (i'i + \v Z u 0 )kl - av 2 ko ~ ~ ^^kl ] 2 - 

a 3 ^i(^-^ 2 )] 2 + [a 3 { Vl -\k^)p 2 h] 2 + €>{a 7 ). (B40) 

Substituting (B 1) and (B2) into (B39) and the result into (B40) and equating like powers 
of <7 shows that 

(^rO +4 £ ) — X B C r o (B 41 a) 

and 

fcosflH ~ \ &io — ± Im \(^X 2 B cl 0 wi cos29) 2 + (fihsinO) 2 , (B416) 

V cos 6/ l 
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where 9 = arctan(/ 9 / 2 a r0 ) is the obliqueness angle of the instability wave and it follows from 
(B9), (B 14) and (B31) that 


<X>1 


1 y+^Ty 00 u . - 

nX%J- n [Jo l(g-9c) 3 9 + 1 " C _ 


dC- 


(B 42) 


Setting h — 0 shows that, apart from differences in notation, (B41) reproduces the results 
given in GW, i.e. their (5.30) and (5.34). In the present analysis, the ± in (B416) correspond 
to solutions for p that are either odd or even functions of C depending on the particular choice 
for U. 

Away from 7 = i (B 40) reduces to 


d—\dQ(i\ + 0 (< 7 ^) — 0 . 


(B 43) 


The solutions to (B43) that match onto (B41) as 7 


x are 


K + Ar-D 2 ]‘ /2 = ^>, (cos^ + = U)a!i' = ^_, (B44) 


and 


nr {II)4 

-(ii) _ __£h) 

i0 4A b ' 


(B45) 


(“^^V) ,/! = A»c-> (cosS'") + -U) 

By including the next-order corrections to (B44) and (B45) and then introducing (B35), 
one can show that, when the minus sign is taken in (B416), (B41) matches onto (B44) as 
h — >• -00 and onto (B 45 ) as h — > +00 or vice versa when the plus sign is taken in (B416). 
It follows that (B41), (B44) and (B45) can be combined to form a set of (multiplicative) 
composite equations that are uniformly valid for 0 ^ 7 ^ 5 . These equations read 


“rt > 2 + (? - l) 2 ] 1/2 - A*4o , 


(B 46a) 
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cos 


COS 


( 7)4 


0 <'> + -K-) & in - 


*° 4A* 


Im{[|a 3 A^ o u; 1 cos 2^ 2 + [m-i)sin0] 2 } 1/2 , (B466) 


and 


cos 


0 (//> + 




(d|"' 2 +4V) 

1 ' 

\ 7TC (//)4 

\ :Un _ nC rO 

COS dWj 

) «tO 

4A b 

+ 4( 


| Im {[\c z \ 2 B c\ 

a 3 V 

£r 0 j 



(B47a) 


where c r0 and u>i are determined from (B 41a) and (B 42), respectively. The results predicted 
by (B46) and (B47) are easily extended to the range | ^ 7 ^ 1 by using the relation (5.4). 
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